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ABSTRACT 

Atmospheric  correction  of  IR  satellite  data  for  retrieval  of 
sea  surface  temperatures  (SST)  by  means  of  differential  absorp¬ 
tion  can  take  on  two  forms.  The  most  popular  method  utilizes 
changes  in  the  spectral  transmittance  of  the  atmosphere  between 
two  or  more  IR  channels  (multichannel).  The  second  form  employs 
path  length  differences  in  the  same  spectral  band  using  two  sat¬ 
ellites.  The  multichannel  and  two  satellite  approaches  are  ap¬ 
plied  to  coincident  imagery  and  in  situ  data  and  compared.  In 
this  test  case  the  multichannel  method  resulted  in  a  mean  SST 
error  of  -3.3°C  vrfiile  the  two  satellite  algorithm  produced  a  mean 
error  of  -0.4°C.  Based  on  radiative  transfer  calculations  with 
the  L0WTRAN-5A  model,  it  is  hypothesized  that  atmospheric  aero¬ 
sols  are  responsible  for  the  poor  performance  of  the  multichannel 
method. 
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I  LUSTRATIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Table  1. 


The  percent  transmittance  of  radiant  energy  as  a  3 
function  of  wavelength  for  the  spectral  range 
which  includes  the  IR  channels  on  the  present 
NOAA  series  of  meteorological  satellites. 

Atmospheric  attenuation  at  AVHRR  channel  4  wave-  i 
lengths  (11  vs.  attenuation  at  channel  3 
wavelengths  (3.7  pm)  based  on  L0WTRAN-5A  calcu¬ 
lations  for  various  model  atmospheres  with  and 
without  maritime  aerosols.  Solid  curve  is  best 
fit  resulting  from  similar  work  by  Deschamps  and 
Phulpin  (1980).  Region  to  the  left  of  the  dotted 
line  is  the  rejection  region  for  the  NOAA-NESS 
Low  Stratus  Test. 

Error  in  MCSST  vs.  scan  angle  based  on  L0WTRAN-5A  7 
calculations  for  a  mid-latitude  winter  atmos¬ 
pheric  model  with  a  surface  temperature  of 
22°C.  Errors  are  plotted  for  both  the  NOAA-NESS 
(McClain,  1981)  and  the  Deschamps  and  Phulpin 
(1980)  algorithms. 

L0WTRAN-5A  calculated  spectral  transmittance  of  7 
the  1962  U.S.  Standard  Atmosphere  with  and  with¬ 
out  inclusion  of  the  maritime  haze  model. 


TABLES 

MCSST  and  TSSST  comparisons 
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COMPARISON  OF  MULTICHANNEL  AND  TWO-SATELLITE  METHODS  FOR  REMOTE  MEASUREMENT 

OF  SEA  SURFACE  TEMPERATURE 


I.  INTRODUCTION 


For  years,  the  problen  of  atmospheric  corrections  has  limited  oceanographic 
uses  of  infrared  (IR)  imagery  to  qualitative  identification  of  ocean  thermal  fea¬ 
tures.  The  atmospheric  attenuation  of  the  IR  signal  effectively  suppresses  the  ac¬ 
tual  sea  surface  temperature  (SST)  and  the  SST  gradients.  Over  the  last  decade, 
scientists  have  addressed  the  problem  by  postulating  or  testing  the  use  of  multiple 
IR  channel  data  to  adjust  for  the  intervening  atmosphere  (Anding  and  Kauth,  1970; 
Prabhakara  et  al.,  1974;  McMillin,  1975;  Deschamps  and  Phulpin,  1980;  and  others). 
The  National  Earth  Satellite  Service  (NESS)  has  performed  atmospheric  corrections 
and  produced  operational  quantitative  SST  products  since  1972  (McClain,  et  al . 
1982),  but  the  early  efforts  were  not  completely  satisfactory. 

Only  recently  have  accurate  quantitative  measurements  resulted  from  IR  digital 
imagery.  The  Advanced  Very  High  Resolution  Radiometers  (AVHRR)  on  TIROS-N,  NOAA-6, 
and  N0AA-7  all  were  configured  with  multiple  IR  channels  and  high  thermal  resolu¬ 
tion.  The  different  atmospheric  transmi ttances  in  these  channels  have  enabled  scien¬ 
tists  to  correct  the  IR  signal  that  is  atmospherically  attenuated,  largely  by  water 
vapor  and  aerosols.  Successful  tests  of  multiple  window  atmospheric  corrections  with 
NQAA-7  data  at  the  National  Oceanic  and'  Atmospheric  Administration  (NOAA)-NESS 
(McClain,  1981)  prompted  the  use  of  multichannel  SST  (MCSST)  methods  in  NOAA  SST 
analysis  beginning  November  1981. 

Another  approach  to  th*'  IR  attenuation  problem  is  to  utilize  two  satellites 
with  at  'ast  one  IR  channe1  viewing  the  same  ocean  surface.  This  possibility  was 
examined  theoretically  by  McMillin  (1975).  Holyer  (1982)  has  proven  that  atmospher¬ 
ic  corrections  can  be  applied  when  coincident  Geostationary  Operational  Environmen¬ 
tal  Satellite  (GOES)  and  NOAA-6  images  are  available.  This  report  summarizes  both 
theories,  and  reports  the  results  of  the  first  quantitative  comparison  between  the 
methods.  An  attempt  is  then  made  to  explain  the  unexpected  results  of  this  compari¬ 
son. 


1 1 .  VERIFICATION  DATA 


The  comparison  data  used  for  this  study  incorporates  multichannel  IR  data  from 
NOAA-6  (Channels  3  and  4)  and  10.5-12.5  pm  data  from  GOES-East  on  February  4,  1980, 
over  the  Gulf  of  Mexico.  The  polar  orbiter  pass  at  1308Z  (7:08  AM  CST)  revealed  the 
northern  Gulf  had  clear  skies,  while  partly  cloudy  conditions  existed  southward  over 
the  Loop  Current.  A  corresponding  G0ES-E  IR  image  at  1200Z  was  processed  for  this 
analysis. 

Surface  truth  for  the  comparison  was  obtained  from  a  variety  of  platforms.  A 
NAVOCEANO  P-3  flight  provided  accurate  SST  and  subsurface  thermal  structure  data  via 
airborne  expendable  bathythermographs  (AXBT)  and  a  Precision  Radiation  Thermometer 
(PRT-5).  NOAA  Data  Buoy  Office  (NDBO)  buoys  and  National  Marine  Fisheries  Service 
(NMFS)  COOP  XBT  and  Japanese  tuna  boat  bucket  temperatures  increased  the  available 
verification  SST  observations  to  64.  However,  cloud  screening  procedures  eliminated 
many  of  these  possible  surface  truth  points  that  occurred  within  12  hours  of  the 
satellite  overpass. 

Obvious  cloud  or  land  contaminated  picture  elenents  (pixels)  were  deleted  by 
first  restricting  the  usable  range  to  10-30°C.  This  range  encompasses  all  valid 


Gulf  of  Mexico  SSTs  except  for  those  in  extremely  shallow  water,  especially  near  the 
mouth  of  the  Mississipi  River.  At  those  surface  truth  locations  where  the  satellite 
temperatures  passed  the  range  test,  a  5x5  array  of  pixels  centered  on  the  surface 
truth  location  was  then  screened  in  NOAA-6  channel  4,  so  that  only  those  that  con¬ 
tained  a  variance  of  less  than  0.3°C  about  the  mean  were  retained.  Only  seven  of  the 
possible  64  surface  truth  points  survived  these  screening  procedures. 

The  dual-window  MCSST  algorithms  evaluated  here  are  normally  for  nighttime  data 
only,  since  AVHRR  channel  3  is  susceptible  to  sun  glint.  The  nadir  track  of  orbit 
3154  was  east  of  the  test  area,  thus  the  sensor  was  scanning  west  toward  the  area  of 
interest.  This  westward-looki ng  orientation  should  avoid  all  possibility  of  sun 
glint  contamination  even  though  this  early  morning  pass  was  in  daylight  conditions. 
Calculations  using  the  visible  NOAA-6  channels  indicated  the  reflectance  values  at 
the  locations  of  the  seven  remaining  surface  truth  points  were  <5.0%  in  each  case, 
thus  confirming  the  absence  of  sun  glint.  These  "cloud  free"  surface  points  serve  as 
verification  points  for  testing  MCSST  and  two  satellite  SST  (TSSST)  algorithms. 

It  should  be  noted  as  background  information  for  this  accuracy  evaluation  study 
that  the  calibration  accuracy  of  NOAA-6  has  been  in  doubt  since  Otis  Brown  at  the 
University  of  Miami  reported  calibration  problems  shortly  after  launch.  Subsequent¬ 
ly,  NOAA-NESS  issued  the  following  statement  in  February  1981  in  APT  81-1: 

Those  individuals/organizations  concerned  with  the  precise  thermal  cali¬ 
brations  should  be  aware  that  a  2  degree  K  thermal  gradient  exists  across 
the  backscan  housing  area  of  the  NOAA-6  AVHRR  instrument.  (The  backscan 
housing  area  is  used  by  the  AVHRR  as  a  blackbody  target  for  thermal  cali¬ 
bration).  When  calculating  thermal  calibration  coefficients  a  simple 
arithmetic  mean  of  the  four  PRT  sensors  is  representati ve  of  the  target 
area  temperature  and  can  be  used  in  these  calculations. 

III.  MULTICHANNEL  METHOD 

The  principle  behind  the  multichannel  correction  is  illustrated  in  Figure  1. 
The  two  IR  channels  aboard  NOAA-6  (3.55-3.93  m  m)  and  (10.5-11 .5  m  m) ,  referred  to  as 
3.7  Mm  and  11  m  m,  respectively,  fall  within  "atmospheric  windows."  These  spectral 
bands  thus  have  relatively  high  transmittances  with  regard  to  middle  and  far  infra¬ 
red  energy  emitted  by  the  ocean  surface.  The  most  significant  atmospheric  absorption 
constituents  in  these  regions  are  water  vapor  and  aerosols.  However,  the  3.7  m  m 
channel  has  a  transmittance  of  ~90%  compared  to  ~75%  for  the  11m  m  channel .  This 
variation  in  infrared  transmittance  between  the  two  channels  (differential  absorp¬ 
tion)  allows  a  correction  to  be  made  for  atmospheric  effects  (Anding  and  Kauth, 
1970). 

High  noise  levels  in  the  TIROS-N  channel  3  data  precluded  any  real  use  of 
multichannel  algorithms  on  this  data.  NOAA-6  multichannel  data  has  been  used  for  re¬ 
search  purposes  (detailed  later)  but  was  never  used  operationally  due  to  software 
limitations.  Offline  tests  with  N0AA-7  data  at  NOAA-NESS  in  1981  proved  that  accu¬ 
rate  SST  retrievals  (RMS  <  1.0°C)  could  be  derived.  These  promising  rerjlts  con¬ 
vinced  NESS  to  operationally  replace  single  IR  channel  derived  SSTs  (corrected  with 
data  from  the  atmospheric  sounder)  in  the  G0SSTC0MP  analysis  with  MCSSTs  on  November 
17,  1981. 

A  large  verification  program  is  now  in  effect  at  NESS  using  fixed  and  drifting 
buoys.  Early  results  show  that  a  global  RMS  of  <1.0°C  is  being  reached  with  a  small 
(-0.5°C)  bias  on  the  cold  side  (McClain,  et  al .  1982).  These  positive  results 
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Figure  1.  The  percent  transmittance  of  radiant  energy  as  a  function  of  wavelength 
for  the  spectral  range  which  includes  the  IR  channels  on  the  present  NOAA  series 
of  meteorological  satellites. 


Figure  2.  Atmospheric  attenuation  at  AVHRR  channel  4  wavelegnths  (11  jim)  vs.  atten¬ 
uation  at  channel  3  wavelengths  (3.7  jim)  based  on  L0WTRAN-5A  calculations  for 
various  model  atmospheres  with  and  without  maritime  aerosols.  Solid  curve  is  best 
fit  resulting  from  similar  work  by  Deschamps  and  Phulpin  (1980).  Region  to  the  left 
of  the  dotted  line  is  the  rejection  region  for  the  NOAA-NESS  Low  Stratus  Test. 
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prompted  NOROA  Code  335  to  encourage  Fleet  Numerical  Oceanography  Center  (FNOC)  to 
again  include  satellite  SSTs  in  their  SST  analysis.  Jim  Cornelius  and  others  at  FNOC 
made  it  possible  for  this  to  become  operational  May  23,  1982,  using  one  8  km  MCSST 
retrieval  received  from  NESS  in  each  2.5  x  2.5  degree  box  (due  to  slow  data  rate  of 
the  communication  line).  The  following  sections  present  the  first  NORDA  effort  to 
compare  the  MCSST  and  TSSST  algorithms  and  explore  any  discrepancies  that  result. 
This  will  directly  benefit  FNOC  products  and  the  satellite-numerical  ocean  model  in¬ 
teraction  project  that  exists  between  NORDA  Codes  335  and  322. 

IV.  TWO-SATELLITE  METHOD 

The  two-satellite  method  utilizes  near  simultaneous  coverage  by  two  satellites 
(normally  one  geostationary  and  one  polar  orbiter,  as  depicted  on  the  front  cover) 
to  make  a  correction  based  on  temperature  differences  resulting  from  different  at¬ 
mospheric  path  lengths  for  the  satellites  as  they  view  a  given  point.  A  simple 
radiative  transfer  model  of  the  atmosphere  is  used  to  translate  geometrical  parame¬ 
ters  and  temperature  differences  into  atmospheric  corrections.  A  complete  descrip¬ 
tion  of  the  TSSST  method  is  in  press  (Holyer,  1982).  The  first  evaluation  of  TSSST 
utilized  data  from  the  Gulf  of  Mexico  in  December  1978.  The  results  were  good,  but 
it  was  not  possible  to  perform  a  comparison  with  MCSST  because  multispectral  satel¬ 
lite  IR  data  was  not  available  at  that  time.  The  present  data  set  is  the  second 
evaluation  of  TSSST  but  the  first  MCSST/TSSST  comparison. 

V.  COMPARISON  OF  MCSST  AND  TSSST  RESULTS 

The  results  of  the  application  of  the  two  methods  to  the  NOAA-6  data  previously 
described  are  shown  in  Table  1  for  the  seven  data  points  which  passed  the  cloud  con¬ 
tamination  screening. 


Table  1.  MCSST  and  TSSST  comparisons 


Ground  Truth 

■Hf 

MCSST 

AXBT  Temp 

SST  Error 

20.4 

19.9  -0.5 

17.3  -3.1 

21.0 

20.3  -0.7 

17.8  -3.2 

21.3 

20.8  -0.5 

18.0  -3.3 

23.7 

21.0  -2.7 

19.4  -4.3 

24.3 

21.7  -2.6 

20.8  -3.5 

24.4 

22.7  -1.7 

21.6  -2.8 

24.1 

23.0  -1.1 

21.1  -3.0 

Mean  Error 

-1.4 

-3.3 

RMS  Error 

0.88 

0.43 

Expected  Bias  (see  below) 

-1.0 

0.0 

True  Error 

-0.4 

-3.3 

The  NAVOCEANO  P-3  aircraft  that  dropped  the  AXBT s  was  equipped  with  a  PRT-5 
radiometer  to  measure  the  radiometric  skin  temperature  of  the  ocean.  Comparison  of 
radiometric  temperatures  (after  careful  calibration  and  atmospheric  correction)  with 
AX8T  temperatures  shows  the  skin/subsurface  differential  to  average  -1.0°C  for  the 
seven  ground  truth  locations  given  in  Table  1.  The  TSSST  method  retrieves  radiomet¬ 
ric  temperatures,  so  a  mean  error  of  -1.0°C  is  expected  when  TSSST  values  are  com¬ 
pared  with  AX3T  data.  The  MCSST  method,  on  the  other  hand,  is  reportedly  bias  cor¬ 
rected  to  give  a  zero  mean  relative  to  subsurface  contact  temperatures.  Thus,  the 
mean  errors  for  the  TSSST  and  MCSST  methods  are  actually  -0.4  and  -3.3°C,  respec¬ 
tively. 

VI.  DISCUSSION 


The  small  number  of  data  points  precludes  any  precise  quantitative  comparison 
of  TSSST  and  MCSST.  However,  in  this  first  head-to-head  competi tion ,  the  TSSST  algo¬ 
rithm  significantly  outperformed  the  MCSST.  The  -3.3°C  bias  in  MCSST  comes  as  a  sur¬ 
prise  in  light  of  the  fact  that  the  latest  figures  from  NOAA-NESS  (McClain,  et  al . 
1982)  show  a  bias  of  <-0.42°C  based  on  267  comparisons  of  N0AA-7  data  with  drifting 
buoy  observations.  Early  NESS  work  with  N0AA-6  data  indicated  a  bias  of  -0.6°C 
(McClain,  1980). 

As  for  other  satellite  sea  surface  temperature  researchers,  the  University  of 
Miami/Rosenstiel  School  of  Marine  and  Atmospheric  Science  (RSMAS)  is  using  the  same 
NOAA-NESS  algorithm  as  is  used  here,  typically  finding  a  bias  of  -1.0°C  with  NOAA-6 
data  (0.  B.  Brown,  personal  communication,  1982).  Dr.  Brown  also  indicates  that  Dr. 
Bernstein  at  the  Scripps  Institution  of  Oceanography  (SIO),  using  the  dual  window 
algorithms  of  Deschamps  and  Phulpin  (1980),  is  finding  a  bias  of  -0.8°C  for  NOAA-6. 
If  NOAA-NESS,  RSMAS,  and  SIO  are  finding  biases  ranging  from  -0.60  to  -1.0°C,  then 
why  does  the  present  data  result  in  a  -3.3°C  bias? 

The  obvious  first  guess  at  an  answer  to  this  question  is  that  some  error  exists 
in  the  NOAA-6  software  or  in  the  way  the  algorithm  is  being  applied.  However,  care¬ 
ful  checking  of  NORDA  software  has  not  yet  revealed  any  such  bugs.  As  indicated  pre¬ 
viously,  the  calibration  of  the  NOAA-6  AVHRR  instrument  has  been  questioned  in  the 
past.  However,  since  RSMAS  and  SIO  get  reasonable  results  from  NOAA-6  and  TSSST  re¬ 
sults  reported  here  are  fairly  accurate,  it  does  not  appear  at  this  time  that  the 
large  MCSST  bias  can  be  attributed  to  sensor  mi  seal ibration.  This,  then,  rules  out 
all  the  easy  answers. 

VII.  RAO  I AT  I VE  TRANSFER  CALCULATIONS 


The  L0WTRAN-5A  computer  code  (Kneizys  et  al.,  1980)  was  used  to  perform  atmos¬ 
pheric  radiative  transfer  calculations  to  shed  light  on  the  large  negative  MCSST 
bias  problem.  Since  the  Oeschamps  and  Phulpin  (1980),  henceforth  referred  to  as  DP, 
algorithms  were  derived  using  the  LOWTRAN  code  (version  3B  rather  than  5A),  the 
first  step  was  to  try  to  duplicate  the  DP  results.  The  DP  work  involved  the  six 
model  atmospheres  provided  by  the  LOWTRAN  code:  tropical ,  mid-latitude  summer,  mid¬ 
latitude  winter,  sub-arctic  winter,  sub-arctic  summer,  and  U.S.  standard.  For  each 
model  three  surface  temperatures  and  two  scan  angles  were  used  giving  a  total  of  36 
possible  combinations.  Only  17  of  the  possible  36  points  are  reproduced  here,  but 
these  are  enough  points  to  determine  if  the  DP  results  could  be  duplicated.  A  figure 
similar  to  one  from  DP  is  shown  here  as  Figure  2  where  atmospheric  attenuation  in 
AVHRR  channel  3  is  plotted  against  attenuation  in  AVHRR  channel  4.  The  solid  line  in 
Figure  2  is  the  best  fit  to  the  36  points  of  DP.  LOWTRAN  provides  haze  or  aerosol 
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models  as  an  option  in  the  program.  It  was  found  that  in  order  for  L0WTRAN-5A  to  re¬ 
produce  the  OP  line  all  aerosol  calculations  had  to  be  eliminated.  Without  aerosol 
effects,  our  results  for  the  17  atmospheric  combinations  (shown  as  circles  in  Figure 
2)  are  quite  close,  but  not  identical  to  the  line  in  Figure  2.  The  slight  differ¬ 
ence  probably  results  from  the  fact  that  DP  used  "nominal"  spectral  response  curves 
for  the  AVHRR,  while  this  study  used  actual  NOAA-6  values.  Also  some  minor  adjust¬ 
ments  have  been  made  to  LOWTRAN  between  versions  3B  and  5A.  The  present  calculations 
were  close  enough  to  DP  to  conclude  that  (a)  L0WTRAN-5A,  recently  operational  at 
NORDA,  was  working  correctly,  and  (b)  the  DP  algorithms  were  based  on  radiative 
transfer  calculations  that  did  not  include  aerosol  effects. 

The  second  step  in  radiative  transfer  calculations  was  to  include  the  effects 
of  aerosols  by  invoking  the  maritime  haze  model  in  L0WTRAN-5A  and  repeating  the  17 
points  in  Figure  2.  The  maritime  haze  model  results  are  shown  with  x's  in  Figure  2 
where  they  differ  significantly  from  the  no-haze  case.  Could  aerosols  in  the  Febru¬ 
ary  4,  1980,  data  be  responsible  for  the  -3.3°C  bias? 

The  L0WTRAN-5A  code  was  used  to  simulate  the  February  4  data  by  using  the  mid¬ 
latitude  winter  atmospheric  model  and  specifying  the  surface  temperature  to  be  22°C. 
The  apparent  surface  temperatures,  as  sensed  from  satellite  altitude  in  AVHRR  chan¬ 
nels  3  and  4,  were  calculated  for  this  situation  with  and  without  maritime  haze  for 
scan  angles  from  0  to  50  degrees.  These  simulated  tanperatures  were  input  to  both 
the  NOAA-NESS  and  DP  MCSST  algorithms  to  retrieve  an  atmospherically  corrected  sea 
surface  temperature.  The  results  of  this  application  of  MCSST  algorithms  to  simu¬ 
lated  NOAA-6  data  are  shown  in  Figure  3  where  MCSST  error  is  plotted  as  a  function 
of  scan  angle  for  no  haze  and  maritime  haze  cases.  From  Figure  3  the  following  con¬ 
clusions  can  be  drawn: 

•  Both  NOAA-NESS  and  DP  algorithms  give  similar  results. 

•  Without  aerosols,  the  MCSST  algorithms  give  accurate  results  with  little 
angular  dependency. 

t  With  a  maritime  aerosol  model  included,  both  algorithms  show  a  large 
negative  error. 

•  When  aerosols  are  included  a  strong  angular  dependency  is  exhibited  as  the 
expected  bias  ranges  from  -2.0°C  at  nadir  to  greater  than  -6.0°C  near  the 
ends  of  the  scan  line. 

•  Aerosol  effects  can  be  shown,  via  radiative  transfer  calculations,  to  pro¬ 
duce  MCSST  errors  comparable  to  the  observed  bias  in  the  February  4,  1980, 
NOAA-6  data  which  did,  incidently,  come  near  the  end  of  the  AVHRR  scan  line. 

VIII.  CHARACTERISTICS  OF  THE  MARITIME  AEROSOL 


A  complete  discussion  of  the  properties  of  marine  aerosols  is  beyond  the  scope 
of  this  report,  but  one  aerosol  property  should  be  discussed.  Figure  4  shows  the 
transmittance  (surface  to  space  at  a  45°  angle  with  respect  to  vertical)  of  the  U.S. 
Standard  Atmosphere  with  and  without  a  maritime  haze.  Note  that  the  transmittance  at 
AVHRR  channel  3  wavelengths  goes  down  relative  to  channel  4  wavelengths  with  the  ad¬ 
dition  of  the  aerosols.  This  behavior  is  opposite  the  spectral  behavior  of  water 
vapor,  which  attenuates  more  in  channel  4  than  in  channel  3.  In  other  words,  the 
maritime  aerosols,  if  present,  will  tend  to  cancel  out  the  effect  of  the  water 
vapor;  hence,  the  dual-window  MCSST  correction  method  has  a  major  drawback. 
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Figure  3.  Error  in  MCSST  vs.  scan  angle  based  on  L0WTRAN-5A  calculations  for  a  mid- 
latitude  winter  atmospheric  model  with  a  surface  temperature  of  22°C.  Errors  are 
plotted  for  both  the  NOAA-NESS  (McClain,  1981)  and  the  Deschamps  and  Phulpin  (1980) 
algorithms. 


Figure  4.  L0WTRAN-5A  calculated  spectral  transmittance  of  the  1962  U.S.  Standard 
Atmosphere  with  and  without  inclusion  of  the  maritime  haze  model. 


Why  then  do  NOAA-NESS,  RSMAS,  and  $10  get  reasonably  good  MCSST  retrievals.  We 
don't  know  about  RSMAS  and  SIO,  but  NOAA-NESS  applies  a  number  of  “cloud  screening" 
tests  to  the  data  (McClain,  et  al.,  1982)  that  also  appear  to  screen  out  high 
aerosol  cases.  Two  NOAA-NESS  tests  seen  especially  relevant  to  the  aerosol  problem. 
First,  the  Low  Stratus  Test  states  that  if  the  temperature  in  channel  4  is  more  than 
0.7°K  higher  than  the  temperature  in  channel  3,  the  data  is  rejected.  This  test  is 
represented  by  the  dotted  line  in  Figure  2,  where  the  area  to  the  left  of  the  dotted 
line  in  Figure  2  is  the  rejection  region.  It  can  be  seen  that  the  LOWTRAN  calculat¬ 
ions  predict  that  some  of  the  clear  sky  with  maritime  aerosol  cases  will  oe 
eliminated  by  the  Low  Stratus  Test.  The  second  important  test  applied  by  NOAA-NESS 
to  all  daytime  data  is  the  Visible  Cloud  Test,  which  states  that  if  the  bidirection¬ 
al  reflectance  measured  by  channel  2  (near-IR  wavelengths)  is  >1.7%,  the  data  is 
rejected.  All  seven  points  analyzed  here  would  have  failed  the  Visible  Cloud  Test 
even  though  they  were  cloud  free.  Thus,  the  NOAA-6  data  processed  here  would  have 
been  rejected  by  the  NOAA-NESS  operational  software  presumably  because  of  cloud 
contamination,  but  probably  because  of  slightly  increased  near-IR  reflectance 
resulting  from  the  presence  of  maritime  aerosols.  It  is,  therefore,  presently 
believed  that  NOAA-NESS  is  maintaining  good  accuracy  in  spite  of  the  aerosol  problem 
by  using  very  restrictive  cloud  tests  that  are  also  eliminating  data  influences  by 
aerosols.  This  procedure  is  necessary  to  get  good  MCSST  results  and  is  acceptable  is 
long  as  the  problem  is  recognized.  Users  must  realize  that  the  price  to  be  paid  for 
the  aerosol  screening  and  resultant  high  accuracy  is  that  one  nay  obtain  a 
relatively  smal 1  number  of  retrievals  even  from  areas  that  are  actually  cloud  free. 

IX.  CONCLUSIONS 


Any  MCSST  algorithm  that  is  based  on  differential  absorption  in  AVMRR  channels 
3  and  4  will  be  adversely  affected  by  the  presence  of  maritime  aerosols.  The  nature 
of  the  degradation  will  be  a  negative  bias,  which  can  easily  be  4-5°C  in  magnitude. 
The  only  way  to  avoid  this  problen  is  to  screen  the  data  to  eliminate  pixels  influ¬ 
enced  by  aerosols  via  some  very  restrictive  "cloud  screening  criteria."  This,  of 
course,  throws  out  an  unknown  number  of  cloud-free  data  points,  thus  reducing  the 
number  of  possible  retrievals.  The  TSSST  method,  since  it  is  based  on  increased  at¬ 
tenuation  with  path  length,  is  not  adversely  affected  by  aerosols.  As  a  matter  of 
fact,  it  does  not  know  the  difference  between  aerosols  and  water  vapor.  The  TSSST 
algorithm  simply  corrects  for  whatever  is  there  as  long  as  the  effects  are  a  func¬ 
tion  of  path  length,  as  described  in  the  TSSST  radiative  transfer  model. 

This  conclusion  is  in  agreement  with  the  recent  work  of  Chedin  et  al .  (19S2) 
where  simulations  of  a  two-satellite  method  were  performed  based  on  METEOSAT  and 
TIROS-N  High  Resolution  Infrared  Sounder  channel  3.  Chedin  et  al .  conclude: 

(The  two  satellite)  way  of  obtaining  different  amounts  of  atmospheric 
absorptions  has  an  important  advantage  over  the  more  usual  way  which  consists 
in  varying  the  wavelength.  This  advantage  canes  from  the  fact  that  the  spectro¬ 
scopic  aspect  of  the  problem  has  been  almost  completely  eliminated.  Indeed,  the 
"anomalous"  (or  unknown)  spectroscopic  behavior  of  some  of  the  components  of 
the  total  absorption  is  still  a  source  of  trouble  in  a  wavelength  differential 
approach. 
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